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Abstract 



pd capture processes at various energies have been analyzed based on solu- 
tions of 3N-Faddeev equations and using modern NN forces. The application 
of the Siegert theorem is compared to the explicit use of ir- and p-like exchange 
currents connected to the AV18 NN interaction. Overall good agreement with 
cross sections and spin observables has been obtained but leaving room for 
improvement in some cases. Feasibility studies for 3NF's consistently included 
in the 3N continuum and the 3N bound state have been performed as well. 
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I. INTRODUCTION 



Recently new nucleon-nucleon (NN) potentials have been worked out, which describe 
the rich NN data base below the pion threshold perfectly well. This set of new potentials, 
often called a new generation of interactions, comprises Nijm I and Nijm II jl], AV18 @ 
and CD-Bonn 0. Used in three-nucleon (3N) scattering calculations they describe many 
observables rather well, especially at lower nucleon lab energies below about 100 MeV Q. An 
important insight thereby is the robustness of that picture: the four NN force predictions 
are very close together. This allows, in case of discrepancies to data, to think that very 
likely one will see three-nucleon force (3NF) effects. Indeed there is increasing evidence that 
certain discrepancies of data and NN force predictions can be cured by adding 3NFs ||. 
Good examples are the minima of the differential cross sections in elastic pd scattering 
at intermediate energies |J or the deuteron vector analyzing power in the same process 0. 
There are also counter examples which clearly demonstrate that the correct spin structure of 
3NFs has not yet been fully established. A most prominent example is the analyzing power 
in elastic pd scattering at low and intermediate energies [f||§. In any case these results 
form a solid basis to study electron scattering and photodisintegration of 3 He or pd-capture 
processes, since the dynamics before and after the photon has been absorbed is fairly well 
under control. Though 3NF effects are of great interest and very likely certain signatures 
have already been identified as pointed out, by far the most 3N observables can to a large 
extent be well described by NN forces only. Also in view of applications to electromagnetic 
processes it is very important to say that the 3 He bound state as well as the 3N scattering 
states can reliably be calculated for one and the same Hamiltonian based on that set of new 
realistic NN forces. 

Now the new dynamical ingredient for electromagnetic processes is the nuclear electro- 
magnetic current operator. Its dependence on the electromagnetic nucleon form factors in 
case of the single nucleon current and its consistency to the chosen nuclear forces are of cen- 
tral interests. Thus the nuclear matrix elements depending on the electromagnetic current 
operator allow the access to the electromagnetic neutron form factors in case of polarized 
3 He targets [|T^-|T2[ and mesonic exchange currents (MECs) provide important insight into 
the nuclear dynamics. 

In this article we concentrate on effects of MECs for the 3N system, which for the bound 



and scattering states is fully treated by solving corresponding Faddeev equations [13,14]. Our 



techniques to handle MECs adapted to our general formalism is described in |T5|]. For low 



photon momenta and processes with real photons the Siegert theorem in a long wavelength 



approximation [16] is quite popular to include some of the two-body currents. We shall also 



use it in a form which does not rely on such an approximation and we shall compare those 
results to calculations using directly MECs. As expected the comparison will look different 
for low and high energies. 

In the long wavelength approximation the Siegert theorem allows one to write the inter- 
action Hamiltonian in terms of the strengths of electric and magnetic fields, i.e., in a gauge 
independent form. The Siegert theorem has been extended before |I7j beyond the long wave- 
length approximation. Diverse aspects and properties of the Siegert-like transformations are 
studied in the above mentioned articles and in 



We shall compare our theory to several recent and older data. Among them are pre- 
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cise analyzing power data taken at Kyushu University tandem accelerator laboratory by 
using an intense <i-beam and a hydrogen gas target of fairly high pressure and by detect- 
ing 3 He particles instead of 7 rays in order to simultaneously measure the whole angular 
distribution |L9| . 

In theory this field has been actively investigated before by many groups. For the older 
calculations we refer to J20| for an overview. More recently Torre |2l] performed for a pd 
capture reaction a very advanced calculation treating initial and final 3N states correctly and 
including MEC's. Ishikawa and Sasakawa [22] even included 3N forces for such a process. 
Further, based on finite rank approximations of NN forces Fonseca and Lehman |?3] studied 
intensively various pd capture reactions. At very low energies there exist highly advanced 
pd capture calculations by Friar et al. [24]] and Kievsky et al. [p5[] , which has been reviewed 



recently in ||26|| . More recently new investigations on photodisintegration and pd-capture 
appeared for very low energies in |27| . [28| . 

Our aim is to use modern forces, keep at least approximate consistency between MEC's 
and the nuclear forces according to the Riska prescription p9| , analyze more data also at 
intermediate energies and compare the results based on Siegert theorem and on the explicit 
use of MEC's. We also do not restrict ourselves to low multipoles like in most of the previous 
work. 

In section [TJ we rephrase the Siegert theorem working only in momentum space. We 
think this is a very transparent and algebraically simple notation in contrast to the rather 
tedious algebraic steps usually presented in a configuration space notation p0| . It also 



does not require long wavelength approximations. In that context we also want to show the 
connection between the partial wave representation of our former work [[JIJ and the multipole 
expansion. There are many kinds of MEC's. In this article we concentrate on ir- and p-like 
exchange currents linked to the specific AV18 NN force. This is outlined in section [ITT]. Then 
in section [TV] we show our way to include a 3NF force into pd capture calculations. The 
broad set of numerical results compared to data is shown in section |V]. In a future work 
we shall compare the different modern NN force predictions with properly related MEC's 
in order to investigate whether this dynamical scenario is robust against interchanges of the 
forces. Finally we summarize in section [VTj. 



II. SIEGERT THEOREM AND MULTIPOLE EXPANSION 



Multipole expansions have a long history and are a natural tool to characterize radiative 
transitions between the numerous levels of nuclei. In the 3N system there is only one bound 
state for 3 He( 3 H) and thus no obvious need for that notation. Also working not only at very 
low energies very many multipoles are involved and we hardly used that notation up to now, 
see however 



32 



The Siegert theorem |16| is embedded in that notation of multipoles and 
that theorem is very powerful in the description of pd capture or the photodisintegration 
processes. Therefore in order to apply it we want to exhibit the connection between our 



partial wave decomposition used up to now [^TJ to the multipole expansion. At the same 
time we want to present a short outline leading to Siegert's theorem carried through in 
momentum space. Work along this line has been presented before in [P3LI33" . 
The nuclear matrix element for photodisintegration of 3 He is 
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JAQ) =< Wf^r^) • 7(o)|*3^ > = ■ T(^), (i) 

where \l/3# e and ^ ^ are 3N bound and scattering states, ~P and P' the total 3N momenta 
before and after photon absorption, ~t%{ Q ) a spherical component of the photon polarization 
vector and j(0) the nuclear current operator. It consists in general of a single nucleon part 
and more than one-nucleon parts. The three components of the nuclear matrix element 
I (Q) can be expanded into spherical harmonics. Choosing the photon direction Q to 
point into the Indirection we expand the product of the polarization vector 1^(1?) and the 
spherical harmonics Yi (Q') into vector spherical harmonics. This then leads immediately to 

t^f) ■ f(\^\z) = -±= £ J2 V2T+1 C(11JM) [dQ'^ m (Q>)-7(\~$\Q>) 



J>\ l=J,J±l 



= -^^^^27+1^(1^1), (2) 

Jcr 

which is the usual multipole expansion. Inserting the Clebsch-Gordan coefficients one gets 
T Ji°{Q) = T%{Q) 

= ~h I d ® { \^7TT^ J+ll( ^ + V ^ ijTi^ JJ - ni ®' )} ' r ^' ) (3) 

and 

TJtiQ) = TjiTiQ) = ^J dQ'^jAQ') ■ T(\-$\Q>) (4) 

This separates the magnetic multipoles Tj^ 19 (/ = J) and the electric multipoles Tf^ (I = 
J ± 1). Now we use the identity |}5f 



QYj,(Q) = "Vlr^imi(O) + V 27Tl^ J - ll( *' (5) 

which allows to rewrite the term with 7 > jj_h(<0) in terms of 7 > jj + n(<0) and most impor- 
tantly the term QYj^(Q). On the other hand Q ■ l(Q) occurs in the continuity equation 
for the electromagnetic current 

V ■ t(V) = < i r * i f ) \[H,p(o)]\*3 He p >=< -p'^KE'm - mm^He'p > 

= U< 'P'¥ f 7 ) \p(0)\^He'? >= QP(V) (6) 

with 

\Q\=Q = u. (7) 
Therefore we end up for the electric multipoles with 

t j\{Q) = 



~ J d&d^j^jj+uifr) ■ + J^Yjs(Q>)p(\-$\Q')}. (8) 
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Another derivation of that formula is presented in the appendix. Note that the first term in 
the curly bracket in Eq. (||) is normally neglected in a long wavelength approximation. The 
last term is now the density matrix element, which is believed to be less affected (at least for 
low momenta) by two-body effects. Therefore a single body operator might be a reasonable 
approximation. In this article we use that assumption and keep also the first term, however 
in a single nucleon current approximation. In the result section we refer by "Siegert" to 
such a choice. Thus the identity given in Eq. (|5|) together with current conservation and 
exact 3N eigenstates of the Hamiltonian H allowed to shift effects of the unknown current 
matrix element I (Q) into p( Q ) and into higher multipoles. Clearly this is only the first 
step towards introducing two-body currents. 

Now let us indicate in the example of the single nucleon density operator the connection 
between our partial wave expansion and the multipole expansion. Working in momentum 
space and in a partial wave decomposition for 3 nucleons we use the basis \pqf3 >, where 
p and q are the magnitudes of Jacobi momenta and (3 a set of discrete quantum numbers 
(orbital and spin angular momenta coupled to the total 3N angular momentum and isospin 
quantum numbers) f3T|j . 

The density matrix element has the structure 

p(Q) = ^ < tfJrW > < pqP\p(Q)\**He > (9) 



For our purpose it is sufficient to regard only < pq(3\p(Q)\^3 He >, which is shown in |3T[ to 
have the form for a single nucleon density operator 

<pq(3(JM)\p(Q || z)\m bound (^M') >= 5 M ,M'hi T ,M T ^\[J4i{-^) S+J 



[/W(t, T, M T )F?(Q) + / T, M T )F?{Q)] £ V<^"WV X'U 

P 1 



E 



X 



5^jfeC(Ai,M;0,0,0)y dxP k {x) 



(2A/ + 1)! ,Ve) A - 

1 <PqP'\y bound > 



q x ' 





Ai 




f / A 




{a 
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k] 


\g V 


A'} 



Y.y/g C(\ 2 ,k,g;0,0,0) 

l 

2 



x {u i I) C(J,g,l;M,0,M'). (10) 



We also indicated now in the bra vector explicitely the total 3N angular momentum J and 
its magnetic quantum number M. Here it is assumed that Q = z. On the other hand the 
part of the electric multipole related to p in Eq. (||) is 



Tf( p = -J dQ'^^-Yj^pdQlQ'). (11) 
In order to perform the integral in Eq. ( |TT]) we need the generalization of 



Eq. ( |10D to an arbitrary Q-direction. This is easily achieved by replacing in Eq. ( 10) 
6M,M>C(J,g,±;M,0,M')by 
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J— Yg^M'-M (Q)C(J, g, -; M, M' - M, M') 

V fi 1 * 
Then the integral can be performed and one ends up with 



Tjf = $<* ( f ) \pqP><pq/3\* e }?> 



where < pqj3\&jf > is given by 



<MP\*j?> 



1 1 



>M T ,M T , 



l±±(-l)tC(J,J,±;M,-t,M') 



JVU-l 



4 v^F V 
I A L 
J V A' 

[&\t, T, M T )Ff (Q) + I^(t, T, M T )F?(Q)} 



J L S 
V \ J 



E $i,i'$s,s'8s,S'$t,t' V A'L' 



E 



A . \ (-A : )!(-A,)! ' ^ 
A2 Ai A 



5] jfe C(Ai, jfe, A; 0, 0, 0) C(A 2 , fc, J; 0, 0, 0) 



dxPk (x) 



A J jfe 

< pqP' Abound > 



(12) 



(13) 



(14) 



As a consequence the variable g in Eq. fllQl) has to be identified with J - the multipole order 
of Eq. (|T4j). A corresponding relation can easily be worked out for the remaining parts in 
Eq. (D and Eq. (§. 



III. MODEL OF MESON EXCHANGE CURRENTS 

The study of meson-exchange currents has also a long history. We follow the scheme 



adopted by the Urbana-Argonne collaboration ||26|| , based on the Riska prescription |29 
Dominant contributions arise from the n- and p-exchanges, to which we restrict ourselves in 
this article. The 7r-current consists of the so called seagull and pion in flight parts, which in 
case of a true pion exchange are displayed in Fig. [TJ The well known expressions are 

f ea9ul \vi,V2) = <[r(l) x r(2)] z F^( Pl )(a(l) . fr)B(2) -v^p 2 )(a(2) ■ p 2 )a(l)\, 
r mic (puP2) = <[r(l) x r (2)],if(fi - p 2 ) [g(l) • ft] [ff (2) • p 2 ] ^ " (15) 

where pi and p 2 are defined in terms of the initial and final momenta as 

-> / -> 
Pi = *i - fci 

p 2 = k 2 - k 2 (16) 
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Further F± is the isovector electromagnetic nucleoli form factor. In case of just a pion 
exchange and dropping the strong form factors v n (p) is given as 



vM = f \ N 2 , (17) 

m% + p z 



where f n NN and m n are the pseudovector irNN coupling constant and the pion mass, re- 
spectively. Similarly the p-current is given as 

J12 = J P ,i + Jp,ii + Jp,iii + 3 P ,iv^ 

Pl ~ F2 

3p,ii(Pi,P2) = ~ i[r(l) x f(2)] z F^(Q 2 ) 

(v p (p 2 ) a(l) x («j(2) x p 2 ) - v p ( Pl ) ff(2) x (ff(l) x ft)) 
3p,m(PnP2) = -i[?(l)xf(2)],W) 

VpiP j~ V & l] ((-(1) >< Pi) ■ (-(2) x pi)) (pi - ft) 
Pl ~~ p 2 

w(pi,p 2 ) - - i[f(i) x mi Fim vpip i~ v & l) 

Pi ~ P2 

(a(2) • (ft x p 2 )(<?(l) x pl) + B{1) • (ft x p 2 )(<?(2) x ft)) . (18) 
Again for just a p-exchange one would have 

\2m N J m 2 p + k 2 
v S p (k) = fan-J^, (19) 

where g P NN, K and m p are the vector, tensor pNN coupling constants and the p meson mass, 
respectively. Now according to Riska's prescription one identifies it- and p-like parts in the 
given NN potential and using the continuity equation finds a linkage of the potential to the 
exchange currents. The procedure is as follows 

v«(k) — v PS (k) = 2v tT (k) - v aT (k), (20) 
v p (k) — ► w(k) = v tT {k) + v aT (k) (21) 

and 

=« T (*), (22) 
where the functions vps(k), vy{k), Vy{k) are evaluated as 

= & J Q drr 2 [j (kr) - l]^(r), (23) 

477- roo 

v tT (k) = — / drr 2 j 2 (kr)v tT (r), (24) 
k Jo 



v T (k) = 4tt / drr 2 j (kr)v T (r). (25) 
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The functions v aT (r), v tT (r), v T (r) are taken from the AV18 potential and are the radial 
functions accompanying the spin-isospin, tensor-isospin and isospin operators, respectively. 
It is interesting to look at the potential dependent functions vps(k), ty(fc) and Vy(k) for 
different NN potentials. In Figs. |2|4j| we compare them for the AV14 |36| and AV18 M. 



IV. INCLUSION OF A THREE-NUCLEON FORCE 

In this section we would like to demonstrate our way of including a three-nucleon force in 
pd capture calculations. For the 3N bound state it has been done in jT3| and also fully con- 



verged calculations are available for elastic and inelastic nucleon-deuteron scattering [§]||[7|]. 

We apply the method proposed in |J7J and start directly with the nuclear matrix element 
for the pd capture process: 

JHQ) =< * 3He l*|-^(3) • ~?(0)|l^ +) > (26) 

The scattering state is then split into three Faddeev components, which for a system 
of identical particles reads 

¥ f +) = fa + fa + fa = (1 + PnP23 + ^23^1 = (1 + P)fa (27) 

The Faddeev component ip fullfils the following equation (we drop the index 1) 

i> = <p + G tPij + (1 + G Q t)G vP{l + P)tJ). (28) 

Here Go is the free 3N propagator, t the off-shell NN t-operator and is a product of the 
deuteron state and a momentum eigenstate of the spectator nucleon. V4 results from the 
decomposition of a 3NF into three parts V} 1 \ which individually are symmetrical under the 
exchange of particles j and k (j 7^ i and k ^ i). 

V 4 = V} 1] + Vf } + Vf } (29) 



Introducing the amplitude T, which obeys the equation [37] 

f = tP<f> + (l + tG Q )vP{l + P)<f> + tPG f + (1 + tG )^ 4 (1) (l + P)G f, (30) 
we end up with 



Eq. (|3"UD reduces to the form: 



^ = + G o f. (31) 



T = tP(p + tPG T, (32) 
when the Hamiltonian contains only NN interactions (V4 = 0). 
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In such a way we study 3NF effects by choosing just either the solution of Eq. (p0|) or 
Eq. Q32D , inserting it into Eq. fl3i"D and Eq. (^7]) and applying the nuclear current operator 
to the proper 3N bound state wave function. 

Thus the nuclear matrix element for the pd capture process S%{Q) can be split into the 
plane wave part and the part containing all initial state interactions 

JtiQ) = < ** He ~P\tt$) ■ 7(0)|(1 + P)0' > 

+ < V*Hj*\t<-{$) • T(0)|(1 + P)G TP*' > (33) 

V. RESULTS 

Let us first repeat calculations performed before by several groups for the magnetic form 
factors of 3 He and 3 H. In Figs. |5|-|6| we compare the single nucleon current prediction with the 
one including the it- and p-like exchange currents. We use the AV18 potential || without 
the various electromagnetic parts. Therefore the result of a Faddeev calculation for the 
3 H binding energy is only 7.623 MeV. The electromagnetic nucleon form factors (Ge and 



Gm) are from ||38|| . We see the well known strong shift of theory towards the data when 
including the MEC's. In this article we restrict ourselves to relatively low photon momenta 
(see Table |), therefore we are not concerned about the discrepancies at the higher Q- values 
in Figs. |5|-|6], where also relativity should play some role. In relation to the higher Q-values 



we refer to p6||39l , where theory has been pushed further including additional MEC's and 
A-admixtures. 

Now we come to the main results and regard first the photon angular distributions for pd 
capture reactions. For the current matrix element J^{Q) defined in Eq. (|I|) the unpolarized 
cross section in the total momentum zero frame (cm) is given as 



da cm = (2vr) V 1 — -1 Yl E \MQ)\ 2 Q 2 dQdn Q 5 (jml He + Q 2 + Q - y/s) 



(34) 



where a, v,ii v p and us are the fine structure constant (~ 1/137), the velocities of the incident 
deuteron and proton, and the energy of the outgoing photon, respectively, y/s is the total 
energy of the p+d and 3 He+7 systems (^/s = Ed + E p = Esne + Q)- The differential cross 
section (da/dQg) is then obtained as 

(d(J r ~ (2vr) 2 «i E E \MQ)? EpEdE3HeQ , (35) 



\dQn J 6 sp v 

where p p is the magnitude of the proton momentum. Note that we use the relativistic phase 
space factor. 

The photon angular distributions for pd capture are shown in Figs. |7j-|l~2"|. The single 



nucleon current prediction is compared to the calculations including the Siegert theorem, 
to the results adding the 7r- and p-like exchange currents to the single nucleon current and 
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to the data. We see that the single nucleoli current prediction underestimates the data. 
Siegert and explicit MEC's are close to the data and to each other, but still leaving a room 
for improvement. At the higher energies when E p = 100, 150 and 200 MeV {E 1 > 70 MeV) 
one may note a slight superiority of the explicit MEC prediction over the Siegert approach. 

In Fig. [13] we show predictions for the 3 He angular distribution in the pd capture process 
for different choices of the electromagnetic current operator at a low energy. Apparently 
the Siegert choice is not fully exhausted by the electric multipole parts of the 7T- and p-like 
MEC's but we see that the 7r-exchange MEC is the dominant one of the two. An interesting 



insight into 3NF effects is shown in Fig. [TJ]. We use the Tuscon- Melbourne force [ £0 , which 
has been adjusted to the 3 He binding energy We see a relatively strong decrease of the 
peak height by including the 3NF in case of Siegert and a less pronounced shift in the same 
direction in case of using MEC's. In |^T| it has been argued (see also j42|) that this fact is 
related to the 3 He binding energy. Apparently this can be more subtle due to the interplay 
with the current operator used. Finally we display the theoretical uncertainty arising from 
the fact that the various NN forces do not yield the same results (Fig. [15]). This spread is 
relatively low (< 5 % in the maximum), which is satisfying in case of Siegert. It will be very 
interesting to see in the future, whether different NN forces taken together with consistently 
applied MEC's will lead to a comparably small spread. Only if such a robustness can be 
demonstrated one can have some confidence in this dynamical scenario. 

Let us now regard various spin observables. Nucleon analyzing powers A y (p) are shown 
in Figs. 16|-|19, At the lowest energy of 5 MeV all three theoretical predictions are close 



together and near the data below 100 degrees, but at extreme forward and backward angles 
the two calculations including two-body currents show strong enhancements, which are not 
seen in the data at the backward angles. Very precise data especially at those extreme angles 
would be of interest to put higher pressure on the theory. At the three higher energies the 
explicit use of MEC's is clearly much closer to the data than the Siegert approach and clearly 
the single nucleon current alone is not acceptable. 

Our result for the deuteron vector analyzing power iT\\ or A y (d) are shown in Figs. p0| 
The agreement with the data is in general only fair. 



22 



At the lowest energy the extreme 
enhancement at large backward angles of the predictions including MEC or Siegert is clearly 
ruled out by the one data point. At E^= 17.5 MeV the MEC prediction appears to be 
reasonable, if there would be not that enhancement at the large angles. At Ed= 95 MeV 
Siegert and the explicit MEC are comparable, but miss the data at backward angles. 

Last not least we regard a rich set of tensor analyzing powers in Figs. |2"3Vf2"9"|. At Ed= 10 
MeV Siegert and MEC predictions agree with each other for T 2 q and T 2 \ and the data. For 
the T 2 2 this is different and the data scatter a lot. We would like to point to the very different 
behavior of the MEC and Siegert predictions at large angles for T 2 q. The reason for that 
unacceptable behavior is right now open. At Ed= 19.8 MeV this quite different behavior for 
MEC and Siegert appears at extreme forward angles. Otherwise both predictions agree with 
each other and the data. At E d = 29.2 MeV our results (0.0326 for Siegert and 0.0315 for 

agree with the calculations by Torre [EH 



MEC) for the single A yy data point (0= 96 °) of 
(0.0339), although in his case the Reid soft-core potential [|3| was used. A yy is fairly well 
described for Ed= 45 MeV and E d = 95 MeV. The data at Ed= 45 MeV have been analyzed 
by us before using the Bonn B potential together with Siegert []32| . 

At E d = 17.5 MeV we show again in detail the various predictions. From Figs. |3"0"|-P^ one 
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can infer that the 7r-exchange provides the strongest shift in relation to the single nucleon 
current prediction, but the additional p-exchange is needed to bring the theory into the 
data. Also Siegert gives essentially the same quality of agreement. Remarkable are again the 
strong enhancements at extreme forward and backward angles which unfortunately cannot 
be checked by available data. 

We would like to add the remark that even at Ed= 17.5 MeV all states for two-nucleon 
total angular momenta up to at least j=2 have to be included. All our calculations are based 
on up to j=3 contributions, which at the higher energies might be not fully converged. In 
Figs. |3"B"|-Po] we display the effect of adding the Tuscon-Melbourne 3NF (adjusted to the 3 He 
binding energy). For all observables the effect is negligible if taken together with Siegert. 
In case of MEC's, however, the shifts are quite noticeable and move the theory somewhat 
away from the data in case of A yy and A zz . 

Finally we demonstrate in Figs. p7|-|40| that the different NN forces taken together with 
Siegert give essentially the same predictions, which is a very nice feature of robustness of 
that dynamical picture. 



VI. SUMMARY 

We have analyzed the angular distributions and some polarization observables in the 
proton-deuteron radiative capture at low and intermediate energies. The corresponding 
photon lab energy for the inverse photodisintegration process would range from 10 to 140 
MeV. We compared the explicit use of MEC's to the use of Siegert's theorem. The MEC's 
were restricted to n- and p-like exchanges as derived from AV18 according to Riska's pre- 
scription. The 3N bound and scattering states are rigorous solutions of the adequate Faddeev 
equations. The calculations are practically converged with respect to angular momentum 
states except possibly at the higher energies. At the lower energies the Siegert and MEC 
predictions are rather close together. The 7r-like MEC provides the strongest shift in relation 
to the single nucleon current towards the Siegert's result, but the p-like piece is important 
as another shift. At the higher energies Siegert and MEC's differ in general, which is to be 
expected since the Siegert approach is mainly active for the lower multipoles. 

The agreement with the data is mostly good but there is some room for improvement. 
Definitely new measurements are needed to improve on certain data sets and to put stronger 
constraints on the behavior of theoretical predictions at extreme angles. Needless to say that 
the correct treatment of the initial state interaction is required in all cases we studied. PWIA 
results are very poor (we did not even display them). 

The strong and sometimes opposite enhancements of the MEC and Siegert predictions 
at extreme angles are worthwhile to be mentioned. Precise data would be very useful to 
rule that out or possibly justify it. In any case if the MEC and the Siegert predictions are 
opposite, the reason for that should be clarified. 

We consider it to be an important result that the different new generation NN force 
predictions together with Siegert are very close together, which demonstrates the stability 
of that dynamical picture. It will be important to investigate in future studies, whether this 
is also true, when explicit MEC's are used consistently to the NN forces. Here we restricted 
ourselves to just one case, when working with MEC, the AV18 potential. 
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Finally we mention that the inclusion of 3NF's in the form of the Tuscon- Melbourne 
27r-exchange (adjusted to the 3 He binding energy) has only a minor effect in conjunction 
with Siegert, but a noticeable one (in some cases) if used together with MEC's. Clearly 
consistency requirements will probably play an essential role and are not the subject of this 
article. 
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APPENDIX A: AN ALTERNATIVE APPROACH TO THE MULTIPOLE 

DECOMPOSITION 

Here we would like to present an independent derivation of Eq. (H), which differs sub- 
stantially in technical tools from the one given in section [Il[ 

The nuclear matrix elements I (Q^ can be expanded in a series of vector spherical har- 



monics 



j(^=£4(Q)?n1(Q), (Ai) 

where 

lUQ) = j dQ'T(Qi)Y} a (Q'). (A2) 
We denote Q = \Q \— Q' . Since / may take values I = | J — 1|, J, J + l in Eq. (|AT|) , we have 
l(Q) = E (iSj-i(Q) Yj1-u{Q) + 4j(Q) Yj1i{Q) + I*jj+x(Q) Yj*hn(Q)) ■ (A3) 



We define 



and 



IjtiQ) = Ijj-i(Q) ?j1-u(Q) + 4j + i(Q) y&Jq) (A4) 



r%"{Q) = 4AQ)?j$(Q). (A5) 

In Eqs. (|A4|) and (|A5|) we separate terms of electric and magnetic type containing the vector 
spherical harmonics with parity (— 1) J and (— 1)( J+1 ) respectively 

=y, ( r &®+ fir (A6) 
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Then we use the following relations [44 



QxY} J+ll (Q) 



J 



2J + 1 



(A7) 



QxYUQ 



j 



2J + 1 



Yjj+ii[Q 



(A8) 



QxYfj-u(Q 



(A9) 



QYjAQ 



J 



V 2J+ 1 
to derive the following identities 



1 J + l 
2J+ 1 



Y 



jj+n\Q 



Y^ 

1 j j +u 



(Q) 



J+l 



2J + 



J 



2J + 1 



QxY 



(A10) 



(All) 



1 jj -li 



J 



l , l ±t±QxYf J1 (Q). 



2J+1 Q J t\~"J " y 2J + 
We use the identities ([All ) and ( A12|) in the angular integrals of Eq. 



$j-i{Q) = JdQ'l(Q')Yfj_ u (Q' 



(A12) 
and get 



J dQ' 



2J 



1 jiQ'-fiQ'^Yj^Q') + iJ^^iQ'xfiQ')]-?}^') 



(A13) 



4j+i(Q)= J dQ'l(Q')Y} J+n (Q') 



dQ' 



J + l 



\[Q'- ?{&)) Yjt{Q') + * \Q' x /(<?')] ■ y} jx (q> 



2J + 



(A14) 



Inserting Eqs. (|A13|) and (|A14| ) into Eq. (|A4l) and making use of Eqs. (|A8|) and (|A10|) leads 
to 



(Q) = QYj\{Q) JdQ'[Q'-T(Q')] Yjz(Q') 



[y;}i(q) x Q] f dQ'[Q' x y J S 1 (g / )]/(5 / 
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QYj*JQ)4:TrT*(Q;l(mg) + i[Yj*j\(Q) x QjAnTf^transv) 



nel , 



(A15) 



In such a way we introduce the electric transverse 



Tf^transv) = £-JdQ'[Q'x Y^^Q')] ■ T(Q' 



4ttQ 

and the electric longitudinal multipoles 

1 



T%{Q-lang) 



4vrQ 



dQ'(Q'-I[Q'))YjAQ'). 



(A16) 
(A17) 

The magnetic multipoles are defined accordingly 

Note that the magnetic multipoles are transverse since 

Q ■ Yjji(Q) = °- The electric lon- 
gitudinal multipoles Tf^{Q]long) do not appear in our studies of the processes with real 
photons. With the help of the identity 



IJ + 1 



J 



QYj^Q) 



'2J + 1 



J 



Q Yj,j+i,i[Q 



(A19) 



one can introduce a longitudinal term into the expression for the transverse electric multi- 
poles. The aim of this transformation is to get a contribution, which can be associated via 
the continuity equation 



q-i(q 



H,p(0)} *s He P 
p(0) V me P) = Qp[Q 



(A20) 



with the matrix elements of the charge density p (q) . So, the electric multipoles can be 
cast into the form 



T%{Q) 



dQ\y/JTlY H {Q)p(Q' 

+ v / 27TTr| J+1 , 1 (g') -i(q' 



(A21) 



which is identical with Eq. (H). 

Finally, choosing the direction of the z-axis along the photon momentum we get for the 
transverse components of the nuclear matrix element 



I(Qe z 



2vr £ V2JT1 \1%{Q) + ZT™ a (Q)} , (£ = ±1 



(A22) 



j>i 



which coincides with Eq. 
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TABLES 



E p 


E d 


Q 


MeV 


MeV 


MeV/c 


5.0 


10.0 


8.8 


8.8 


17.5 


11.4 


9.9 


19.8 


12.1 


14.8 


29.6 


15.4 


22.5 


45.0 


20.5 


47.5 


95.0 


37.2 


100.0 


199.9 


72.3 


150.0 


299.9 


105.7 


200.0 


399.8 


139.1 



TABLE I. Proton and deuteron lab energies and corresponding photon momenta for the inverse 
reaction in the lab system for experiments analyzed in this work. 



18 



FIGURES 




N 



k' 



71 



N 



+ 




K 



N N 



FIG. 1. The " pion- in-flight" and "seagull" pionic exchange currents. 
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FIG. 2. Comparison of the form factors v n for AV18 (solid line) and AV14 (dashed line). 
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FIG. 5. The magnetic form factor of 3 He. The solid line includes MEC (see text), while the 



dashed line is based on a single- nucleon current study. Data are from [45] 
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FIG. 6. Same as in Fig. [| for 3 H. Data are from [45] 
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FIG. 7. The photon angular distribution for pd capture at Ed= 19.8 MeV against the cm. 7-p 
scattering angle. The single nucleon current prediction is given by the dashed line, including MEC 



leads to the solid line and the dotted line is due to Siegert. Data are from [47] 
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0[deg] 

FIG. 13. The 3 He angular distribution for pd capture at Ed= 17.5 MeV against the cm. 3 He-d 
scattering angle. The curves describe the single nucleon current (dashed), the single nucleon plus 
the 7T-MEC (dashed-dotted), the single nucleon plus the ir- and p-MEC (solid) and Siegert (dotted) 
predictions. 




FIG. 14. The 3 He angular distribution for pd capture at Ed= 17.5 MeV against the cm. 3 He-d 
scattering angle. The curves describe the Siegert without 3NF (dotted), the Siegert including 3NF 
(dashed-dotted), the single nucleon plus MEC without 3NF (dashed) and the single nucleon plus 
MEC including 3NF (solid) predictions. 
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FIG. 15. The 3 He angular distribution for pd capture at 17.5 MeV against the cm. 3 He-d 
scattering angle. The curves describe the Siegert predictions based on different NN forces: Nijm I 
(dotted), Nijm II (dashed-dotted), CD Bonn (dashed) and AV18 (soild). 
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FIG. 25. The same as in Fig. 23 for the tensor analyzing power T22. Data are from [46 
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FIG. 26. The same as in Fig. || at E d = 19.8 MeV. Data are from @ (circles) and from |51|| 
(squares) . 
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FIG. 27. The tensor analyzing power A yy at Ed= 29.2 MeV against the cm. 7-d scattering 
angle. Curves as in Fig. pi. The data point is from |21]. 
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FIG. 28. The tensor analyzing power A yy at Ed= 45 MeV against the cm. 7-d scattering 
angle. Curves as in Fig. ^. Data are from J32| . 
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FIG. 33. The tensor analyzing power A xx at -E^ = 17.5 MeV against the cm. 3 He-d scattering 
angle. Curves as in Fig. [14|. Data are from IS]. 



-0.02 - 



-0.1 



— i ^~r~_ — i r 



0.04 - / 



-0.06 - 



-0.08 - 




J I I l L± 



20 40 60 80 r 100 120 140 160 180 

9 [degf 



FIG. 34. The same as in Fig. 33 for the deuteron vector analyzing power A y (d). Data are 



from [191. 
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FIG. 37. The tensor analyzing power A xx at -E^ = 17.5 MeV against the cm. 3 He-d scattering 
angle. Curves as in Fig. Data are from IS]. 
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FIG. 38. The same as in Fig. 37 for the deuteron vector analyzing power A y (d). Data are 



from [191. 
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